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What determines the fraction of elliptical galaxies in 
clusters? 
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ABSTRACT 

We study the correlation between the morphological mix of cluster galaxies and the 
assembly history of the parent cluster by taking advantage of two independently de- 
veloped semi-analytic models for galaxy formation and evolution. In our models, both 
the number of cluster members and that of elliptical members increase as a func- 
tion of cluster mass, in such a way that the resulting elliptical fractions are approx- 
imately independent of cluster mass. The population of cluster ellipticals exhibit a 
marked bimodal distribution as a function of galaxy stellar mass, with a dip at masses 
~ 10 10 M Q . In the framework of our models, this bimodality originates from the com- 
bination of a strongly decreasing number of galaxies with increasing stellar mass, and 
a correspondingly increasing probability of experiencing major mergers. We show that 
the correlation between the measured elliptical fraction and the assembly history of 
the parent cluster is weak, and that it becomes stronger in models that adopt longer 
galaxy merger times. We argue that this results from the combined effect of a decreas- 
ing bulge production due to a reduced number of mergers, and an increasing survival 
probability of pre-existing ellipticals, with the latter process being more important 
than the former. 

Key words: galaxies: formation - galaxies: evolution - galaxies: bulges - galaxies: 
interactions - galaxies: clusters: general 



1 INTRODUCTION 

It has long been known that early type galaxies (ellipticals 
and lenticulars) reside preferentially in dense regions of the 
Universe such as rich clusters, while late type galaxies rep- 
resent a larger fraction of the galaxy population inhabiting 
regions of 'average' density. Such a morphology-density rela- 
' tion was noticed in early observational studies (indications 
of a correlation between the type of nebulae and t he envi- 
ronment can be found in 'The Realm of Nebulae ^bv lHubbld 
1936), and was firmly established by iDresslerl (|l98d ). 

In the past decades, much observational information has 
been collected on the morphological distributions of cosmic 
galax y populations, and o n its depend ence on the environ- 
ment. iButcher fc Oemlerl (|l978al . ll984l ) showed, for the first 
time, that the fraction of blue (star forming) galaxies in 
clusters increases with increasing redshift. Detailed morpho- 
logical studies have been carried out in the following years, 
demonstrating that the fraction of spiral galaxies increases 
with increasing redshift, and that this increase appears ap- 
proximately balanced by a decrease in the fraction of the 
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lenticular galaxies since z ~ 0.5. Over the same redshift 
range, the fraction of elliptical galaxies is appro ximately con- 
stant |Dressler et al. Hl997l : lFasano et al.ll20ob]) . In more re- 
cent years, detailed mo rphological studies h ave been pushed 
to lower mass ranges (| Wilman et al J 120091) . and to higher 
redshift, where the mean fraction of differen t morphological 
types does not appear t o evolve significantly (Po stman et al.l 



types docs not appear t< 
l2005l : lDesai et alj|2007l ) 



At a given redshift, clusters with similar mass (mea- 
sured from either X-ray luminosity, or velocity dispersion) 
exhibit a non neglig ible scatter in t heir morphological com- 
position (e.g. iPoggianti et al.ll2009l ). In the context of the 
currently accepted paradigm for structure formation (the 
ACDM model), it is logical to relate this cluster-to-cluster 
variance to the dynamical history of the cluster. Although 
difficult to test quantitatively, this expectation is confirmed 
by early observations that centrally concentrated clusters 
have typically large populations of ellipticals and lenticulars 
and relatively low numbers of spirals, while irregular, unre- 
laxed clusters are more spiral-rich and show weaker radial 
gradie nts in their morphological mix (e.g. lButcher fc Oemlerl 
Il978bl ). In this paper, we will address this issue by consider- 
ing two different semi-analytic models of galaxy formation, 
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and by relating the predicted fraction of elliptical galaxies 
to the accretion history of the simulated cluster haloes. 



2 THE GALAXY FORMATION MODELS 

In this paper, we take advantage of two independently de- 
veloped galaxy formation model s: the 'Munich' mode l , with 
the implementation discussed in lDe Lucia fc BlaizotI (|2007l ) 
and applied to the Millennium Simulation, and the MOR- 
GANA mode l, as a dapted to the WMAP3 cosmology in 
lLo Faro et ail (|2009T l. Hereafter, we will refer to the former 
model as DLB07. We note t hat in previous wo rk, we have 
used the models presented in Wang et al. (2008) whic h cor- 
responds to the model by |Pe Lucia fc Blaizotj l|2007l ) used 
here, but has been adapted to a WMAP3 cosmology. In this 
paper, we use the model applied to the Millennium Simula- 
tion as this provides a larger volume and therefore a larger 
number of massive haloes. 

The simulations employed in this study assume a dif- 
ferent cosmology: the Millennium Simulation assumed a cos- 
mological model that is consistent with WMAP first-year re- 
sult:^]. As shown in previous work, however, once the model 
is re-tuned to account for the change in cosmolog y, the basic 
resul ts and trends do not change significantly (Wan g et al.l 



2008). In addition, we note that the two models adopt dif- 
ferent definitions for the halo mass: in DLB07, this is given 
by M200 and is computed from the simulation outputs as 
the mass contained in a sphere of radius -R200, for which the 
mean overdensity is 200 times the critical density of the Uni- 
verse at the redshift of interest. For MORGANA, the masses 
are simply given by the sum of the particle mass associated 
with the halo, computed using PINOCCHIO l| Monaco et alj 
120021 ). 

In this section, we provide a brief summary of the model 
elements that are relevant to the present study. We re- 
fer to the original papers for a more detailed discussion of 
the physical processes considered, and of the correspond- 
ing modelling adopted. Both models consider two different 
channels for the formation of bulges: galaxy-galaxy merg- 
ers and disk instabilities. The relative importance of these 
channels, in different envir onments and at differe nt times, 
has been studied in detail in |Pe Lucia et~aH l|201ll ). while in 
iFontanot et al.l (|201ll ) we focused on the statistics and prop- 
erties of bulgeless galaxies. Both models used in this study 
assume a Chabrier Initial Mass Function. 

Mergers are classified as minor or major according to 
their baryonic (gas + stars) mass ratio. If this is smaller 
than 0.3, the merger is classified as minor: the stellar mass 
of the secondary is added to the bulge component of the 
primary galaxy, and the merger is accompanied by a star- 
burst. The resulting stars are added to the bulge component 
(in MORGANA) or to the disk component (in DLB07). If 
the baryonic mass ratio of the merging galaxies is larger 
than 0.3, we assume that we witness a major merger. In 
this case, both models assume that the disk components of 
the merging galaxies are completely destroyed. The remnant 



1 The most important difference between WMAP first and third 
year data is a lower value for the amplitude of matter fluctua- 
tions on 8 h~ 1 Mpc scale (o"g), wh ich leads to a delay in structure 
formation (e.g. IWane et al.ll2008T ). 



spheroidal galaxy can re-grow a new disk, if fed by an ap- 
preciable cooling flow. In previous work, we have found that 
the merger model adopted in MORGANA provides merg- 
ing times that are systematically shorter (by up to an order 
of magnitu de) than those adopte d in the DLB07 model (see 
Section 7 of lDe Lucia et af1l2010r i . The shorter merger times 
adopted in MORGANA lead to a more efficient formation of 
bulges and to larger number densities of early type galaxies, 
particularly at high redshift. As we will show below, the dif- 
ferent modelling adopted for galaxy mergers also affects the 
relation between the morphological fraction and the accre- 
tion history of dark matter haloes. In order to quantify the 
significance of this effect, in the following we will also show 
or discuss the results obtained from MORGANA using the 
same dynamical friction timescale prescription adopted by 
DLB07. 

The treatment of disk instability differs significantly in 
the two models con sidered: both adopt the same stabilit 
criterion proposed in lEfstathiou. Lake fc Negropontel (Il98l 
but use different definitions for the relevant physical quan- 
tities, and make different assumptions about the outcome 
of instabilities: DLB07 only transfer to the bulge a fraction 
of the stellar disc that is enough to restore stability. In the 
MORGANA model, half of the disk baryonic mass (both 
gas and stars) is transfer r ed to the bulge. As discussed and 
shown in lDe Lucia et all l|201ll ). this translates into a more 
relevant contribution of disk instability to bulge formation. 

In the framework of our models, most of the elliptical 
galaxies acquire their morphology through major mergers. 
Disk instability can contribute significantly for low and in- 
termediate mass galaxies, depending on the adopted treat- 
ment for galaxy mergers and instabilities. As mentioned 
above, bulge dominated galaxies can later grow a new disc, if 
they are fed by an appreciable cooling flow. We have shown 
that the rates of disc regrowth are negligible for massive 
galaxies and at low redshift. They represent, however, a non- 
negligible component of the evolution of low and intermedi- 
at e mass galaxies, parti cularly at high redshift (see Section 6 
of lDe Lucia et al.1120111 ). As we focus on galaxy clusters, the 
model ellipticals considered in this paper are almost all satel- 
lite galaxies (with the exclusion of central cluster galaxies). 
For these galaxies, the bulge-to-total ratio is not affected af- 
ter accretion onto a more massive halo in the MORGANA 
model. DLB07 accounts for mergers between satellites (that 
are, however, rare) so that the bulges of satellite galaxies can 
still grow through this physical mechanism. Finally, none of 
the models used in this study include 'environmental' pro- 
cesses such as tidal stripping or harassment, that can poten- 
tially affect the morp hology of satellite galaxie s orbiting in 
a massive cluster (e.g. iMastropietro et al.ll2005T ). 

In our previous work, we have considered alternative 
prescriptions to model bulge formation, including predic- 
tions obtained when the disk instability channel is switched 
off. We have verified that the results presented in the fol- 
lowing do not depend significantly on these assumptions. 
Therefore, we will discuss only results obtained by our de- 
fault models. As these data have not been used to 'tune' 
the models in the first place, they can be considered as gen- 
uine model predictions, and compared with available obser- 
vational measurements. 
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Figure 1. Total numbers of galaxies (left panel), numbers of ellipticals (middle panel), and elliptical fractions (right panel) inside the 
cluster virial radius, as a function of the cluster mass. Filled and open symbols (these are coloured black and red in the online edition of 
the Journal) are obtained considering all galaxies more massive than 10 9 and 10 10 Mq, respectively. 



3 THE FRACTION OF ELLIPTICAL 
GALAXIES IN CLUSTERS 

In this study we have considered galaxy clusters with masses 
in the range logM/M© = 14 — 14.8, at redshift zero. In the 
simulation used by MORGANA, there are 100 clusters over 
this mass range. To compare the predictions from this model 
to those from DLB07, we have selected the same number 
of haloes from the Millennium Simulation, uniformly dis- 
tributed in mass over the same mass range. In the follow- 
ing, we will define as ellipticals all galaxies that have a stel- 
lar bulge to total ratio larger than 0.9. When relevant, we 
will comment on how results can be affected by a different 
threshold. 

Figure [T] shows the total number of galaxies (left panel), 
the number of ellipticals (middle panel), and the elliptical 
fractions (right panels) as a function of the cluster virial 
mass. Only galaxies residing within -R200 have been con- 
sidered in the DLB07 model. Since MORGANA does not 
provide information on the position of galaxies within dark 
matter haloes, we have simply considered in this case all 
galaxies associated with the final cluster. Given the differ- 
ent definitions adopted, any difference between model pre- 
dictions (but we will see that these are very small) should 



be interpreted with caution. Filled and open symbols in Fig- 
ure [1] correspond to galaxies more massive than 10 9 and 
10 10 Mq , respectively. The former limit corresponds to the 
approximate resolution limit of the Millennium Simulation, 
while the latter corresponds to a typical limit for observa- 
tional studies. 

Figure [T] shows that both the total number of galaxies 
and the number of ellipticals increase with increasing halo 
mass. When a lower stellar mass threshold is chosen, the 
predicted numbers are significantly higher. A cluster of mass 
2.5 x 10 14 Mq contains on average ~ 50 (~ 60) galaxies more 
massive than 10 10 M in the DLB07 (MORGANA) model. 
When considering a mass limit of 10 9 Mq , the average num- 
ber of cluster members within the virial radius increases to 
- 190 (~ 260 in MORGANA). The number of ellipticals 
increases too, but not as much as the total number of galax- 
ies. This is expected given that, as the stellar mass increases, 
a larger f raction of galaxies ar e classified as ellipticals (see 
Fig. 7 in |Pe Lucia et al"1l201ll ). Interestingly, the halo oc- 
cupation distribution of the two models used in this study 
is different, with MORGANA always predicting a slightly 
larger number of cluster members with respect to DLB07. 
The difference is significant for the most massive clusters 
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included in our sample, and when considering all galaxies 
more massive than 10 9 Mq. We stress, however, that a dif- 
ferent definition for cluster members has been adopted for 
the two models. In addition, we are using the dynamical 
information available from the simulations to define cluster 
members, while an accurate comparison with observational 
measurements should account for possible contamination by 
interlopers along the line of sight. As mentioned above, the 
merger times adopted in MORGANA are about one order of 
magnitude shorter than those adopted in DLB07. By using 
longer merger times in MORGANA, the number of cluster 
members increases even further. 

The fraction of elliptical galaxies resulting from the 
numbers shown in the left and middle panel of Figure [1] 
does not vary significantly as a function of cluster mass, 
in agreem ent with observational measurements in th e local 
Universe l|Wilman et al.l 120091 ; IPoggianti et alJ liooarl. The 
predicted elliptical fractions are of the order of 10 per cent in 
both models, when all galaxies more massive than 10 9 M 
are considered. For a mass threshold of I0 10 Mq , the ex- 
pected fractions increase, and the scatter becomes larger. 
Interestingly, the halo to halo scatter appears to increase 
slightly with decreasing halo mass. This is more evident 
in the MORGANA model, but we note that in this case 
haloes are not distributed uniformly in mass and the num- 
ber of clusters at the largest masses considered is quite low. 
Therefore, the very narrow range of elliptical fractions pre- 
dicted by this model for the most massive haloes might be 
fortuitous, and just due to poor number statistics. 

When considering all galaxies more massive than 
10 10 M Q , the mean elliptical fraction is 0.22 for the DLB07 
model, and 0.24 for MORGANA. This is lower than the 
average value of ~ 0.32 measured for the Wide-field Nearby 
Galaxy clusters Sur vey (WINGS), using a similar mass cut 
l|Vulcani et al.|[201ll ). We note, however, that only galaxies 
within O.6R200 have been considered in this observational 
study, as this is the largest radius covered in all their cluster 
fields, and that the study is based on a definition of ellipti- 
cals that differs from that adopted in this paper (morpholo- 
gies h ave been assigned u sing V-band images). In previous 
work (tSimard et al.ll2009h . we have shown that the early- 
type fractions predicted by the DLB07 model compare well 
to observational measurements from the Sloan Digital Sky 
Survey (SDSS) in the local Universe and from the ESO Dis- 
tant Cluster Survey (EDisCS) at redshift ~ 0.6. Also in that 
study, however, a different (closer to that used in the obser- 
vations) definition of 'early-type' galaxies was adopted, so 
that the predicted fractions shown in this study are not the 
same as those shown in Simard et al. In a forthcoming pa- 
per (Wilman et al., in preparation), we will carry out a more 
detailed comparison between the observed mix of different 
morphological classes and predictions from our galaxy for- 
mation models. 

The left panels of Figure [2] show the predicted distri- 
butions of stellar masses for all cluster galaxies (thin his- 
tograms) and for the cluster ellipticals (thick histograms). 
These distributions have been obtained by stacking the 



2 We note that both IWilman et al. I j2008h and IPoggianti et al] 
are based on magnitude limited samples, while we are using 
a cut on total stellar mass. 



galaxies in all clusters, and have been normalized to the total 
number of galaxies in each distribution. The two models pro- 
vide very similar predictions, but those from MORGANA 
are more skewed towards less massive galaxies. Interestingly, 
both models predict a bimodal distribution for elliptical 
galaxies, with a pronounced 'dip' around ~ 10 10 Mq. Unfor- 
tunately, this is below or approximately at the limit of the 
observational measu rements for the WINGS sample used in 
IVulcani et all (|201ll ). We note that this bimodal behaviour 
is found in our models also when considering the global 
elliptical population (i.e. not only ellipticals in clusters), 
which does not appear to be supported by available obser - 
vations fe.g. lTrentham fc Hodgkii32002l ; lDriver et alf2003h . 
We stress that our models (like most of the recently pub- 
lished models) overpredict the number densities of small to 
intermediate mass galaxies (|Fontanot et al.l | 200 9) , so the im- 
portance of the peak at small masses is likely over-estimated. 

In both models, the bimodal distribution of the cluster 
elliptical masses is significantly reduced (but still apparent 
in the MORGANA model) when the adopted threshold for 
defining a bulge dominated galaxy as an elliptical is lowered 
to ~ 0.7. In this case, only one peak is visible in the DLB07 
predictions at masses log(M sta r) ~ 10.5. We have verified 
that, in our model, this bimodality is not significantly af- 
fected when the disk instability channel for bulge formation 
is switched off, so a differential efficiency of bulge formation 
through disk instability is not responsible for the shape of 
the cluster elliptical mass distribution sho wn in Figure [2] 
In our previous work (|De Lucia et al.ll201ll ). we have shown 
that disk regrowth is more efficient for intermediate mass 
galaxies. In order to test if this could be responsible for the 
observed dip at intermediate masses, we have calculated the 
mass distribution of all galaxies that have been ellipticals 
in their past, either considering only those surviving at red- 
shift zero (i.e. excluding those that have merged with other 
galaxies) or all galaxies in the merger trees of the cluster el- 
lipticals. In both cases, we find that the predicted mass dis- 
tribution exhibit a marked bimodality, with a pronounced 
'deficit' of elliptical galaxies at intermediate masses. 

We interpret this bimodality as a result of the increas- 
ing probability of sufferi ng a major merger with increasing 
mass (see Figure 9 i n De Lucia et alj|200d and Figure 6 in 
Wang & Kauffmann 2008), and of the strongly decreasing 
number of galaxies of larger masses (as shown by the thin 
lines in Figure[5]). The convolution of these distribution func- 
tions results in a lower number of intermediate mass galaxies 
suffering of major merger events in their past history, com- 
pared to galaxies residing in the low and high mass peaks 
of the distributions sho wn in the left panels of Fi gure [2] 
In particular, Figure 3 of IWang fc Kauffmann! |200cf ) shows 
that the region where the dip in the mass distribution of 
elliptical galaxies is visible, corresponds to a regime where 
the probability of suffering of a minor merger is significantly 
larger than that of experiencing a major merger. This hap- 
pens because, during the time that elapses between a halo 
merger and the actual merger between the galaxies resid- 
ing at the centre of the merging haloes, the stellar mass of 
the satellite does not increase significantly, while the cen- 
tral galaxy grows in mass as it is fed by cooling from the 
surrounding hot halo. As a consequence, the stellar mass 
ratio between the two galaxies decreases, so that a major 
merger between two haloes can lead to a minor merger be- 
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Figure 2. Distribution present day stellar masses (left panel) for model ellipticals in our cluster sample, and for their parent halo mass at 
the time of accretion (right panels). The distributions shown have been obtained by stacking all clusters in the sample, and normalizing 
to the total number in each distribution. Solid and dashed lines refer to the cases when all galaxies more massive than 10 9 and 10 10 M Q 
are considered, respectively. Thin and thick lines (black and red in the online edition of the Journal) are for all galaxies and for those 
that are classified as ellipticals, respectively. 



tween their galaxies. Figure 3 in I Wang fc Kauffmannl ()2008l ) 
shows that the probability of experiencing a minor merger is 
largest at intermediate masses, which explains why lowering 
the adopted bulge-to-total threshold fills the intermediate 
region, and tends to wash out the bimodality. 

For each cluster member in our sample, we have traced 
back their main progenitor until the galaxy is for the last 
time a central galaxy of a dark matter halo, and we have 
recorded the parent halo mass at this time. In the follow- 
ing, we refer to this as the 'time of accretion', although this 
will not always coincide with the time when the galaxy is 
accreted onto the main progenitor of the final cluster (De 
Lucia et al., in preparation). The right panels of Figure [5] 
show the distributions of halo masses at accretion for all 
cluster members (thin lines) and for the ellipticals (thick 
lines). Again, the distributions obtained for all clusters have 
been stacked and normalized to the total number of galaxies 
in each of them. 



The figure shows that, in both models, elliptical galaxies 
tend to be accreted when they reside in more massive haloes 
with respect to the total population of cluster members (the 
distribution predicted by MORGANA has a higher lower 
mass limit than in the DLB07 model: i.e. in MORGANA, 
ellipticals tend to be accreted, on average, in more massive 
haloes than in the DLB07 model). This is not surprising 
considering that ellipticals represent a larger fraction of the 
most massive galaxies, and that there is a relatively tight 
correlation between the galaxy ma ss and that of the parent 
halo mass for central galaxies (e.g. [Wang et al. l l2006l ). 

4 THE ASSEMBLY OF CLUSTERS WITH 
DIFFERENT ELLIPTICAL FRACTIONS 

In the previous section, we have shown that ellipticals tend 
to be accreted in larger haloes, with respect to the entire 
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cluster galaxy population. Given these results, one would 
expect naively that haloes that have acquired a larger frac- 
tion of their mass through the accretion of 'massive' haloes 
would host a larger fraction of ellipticals. One has to con- 
sider, however, that elliptical galaxies can also disappear 
from the sample of cluster members by merging with the 
central galaxies of the hosting halo (or with other satellites 
in the DLB07 model). 

In order to address this issue, we have analysed the 
accretion histories of all clusters included in our sample. For 
each halo, we have traced back in time its main branch, i.e. 
the branch of the tree that is obtained by connecting the 
halo to its mam progenito^. We have then considered all 
substructures residing in the main branch at each time, and 
have traced each of them back in time until they were main 
haloes of a FOF-group. The top panels of Figure[3]show the 
mass distributions of accreted haloes for the clusters that 
host an elliptical fraction lower (thin dashed lines) than the 
10th percentile, and higher (thick solid lines) than the 90th 
percentile of the distribution of elliptical fractions measured 
for all 100 haloes considered. Again, the distributions from 
the two cluster samples have been stacked. Different columns 
correspond to different models, as indicated by the legend, 
while the bottom panels show the corresponding cumulative 
distributions. The differential distributions shown in the top 
panels of Fig. [3] have been weighted by mass in order to 
remove the dominant mass dependency and emphasize the 
differences between the two samples. 

In the DLB07 model (left panels), there is a clear differ- 
ence between the two samples, which is more evident when 
looking at the cumulative distributions: clusters that host 
the highest elliptical fractions also accreted a larger num- 
ber of haloes more massive than ~ 10 11 Mq, with respect 
to the clusters that host the lowest elliptical fractions. In- 
terestingly, the difference between the two samples persist 
over the entire mass range: this implies a lower contribution 
from diffuse accretion for clusters with large elliptical frac- 
tions. In the standard MORGANA model, no significant 
difference is found between the two samples. If, however, 
longer merger times are adopted (right panels), then a dif- 
ference between the haloes with largest and lowest elliptical 
fractions become visible, and it is of the same order of mag- 
nitude of that found in the DLB07 model. At first sight, 
this result appears counter-intuitive because one would ex- 
pect that shorter merger times would translate into a better 
matching between the morphological mix of the galaxy pop- 
ulation and the assembly history of the halo. One has to 
consider, however, that changing the merger times would 
affect the elliptical galaxy population in two distinct ways: 
on the one hand, longer merger times would tend to de- 
crease the number of mergers (and therefore the number 
of bulge dominated galaxies). On the other hand, longer 
merger times would also tend to 'preserve' the pre-existing 

3 In MORGANA, this is simply the most massive progenitor at 
each node of the tree. A different definition is adopted in DLB07: 
the main progenitor is selected by choosing the branch of the 
merger tree that accounts for most of the mass of the final s ystem , 
for the longer period. As explained in lDe Lucia &; Blaizoj l l2007l) , 
this avoids possible problems arising when the selection of the 
most massive progenitor would be ambiguous like, for example, 
when there are two progenitors of similar mass. 



ellipticals from being accreted onto the central galaxies (or 
from merging with other satellites if this physical process 
is included). When adopting longer merger times in MOR- 
GANA, we find that the second process would be slightly 
more important than the first. As a consequence, both the 
total number of cluster members and the number of ellipti- 
cal members would increase. This implies that longer merger 
times preserve a better memory of the accretion history of 
the parent dark matter haloes, thereby creating a stronger 
correlation between the morphological mix of the cluster 
galaxy population and its dynamical status. 



5 DISCUSSION AND CONCLUSIONS 

At fixed cluster mass, the observed properties of the cluster 
galaxy population exhibit a large variation. Such a scatter 
is in part due to observational uncertainties in the observed 
quantities. In the hierarchical framework, however, it is nat- 
ural to link the observed halo-to-halo scatter to a range of 
dynamical histories of the parent cluster. In this paper, we 
have investigated the link between the predicted fraction 
of elliptical galaxies and the accretion history of the parent 
dark matter halo, by taking advantage of two different semi- 
analytic models of galaxy formation. Our main results can 
be summarized as follows: 

• For both models used in our study, the predicted el- 
liptical fractions do not vary significantly as a function 
of the cluster mass, for the range of masses considered 
(M200 ^ 10 14 Mq). This appears to be in qualitative agree- 
ment with observational measurements Jwilman et al.ll2009l ; 
iPoggianti et al.ll2009l : ISimard et al1l2009h . In our models, a 
constant elliptical fraction results from an increasing num- 
ber of both cluster members and elliptical members, as a 
function of cluster mass. 

• Cluster ellipticals exhibit a marked bimodal distribu- 
tion in stellar mass. In both models, the bimodality is 
reduced when a lower (~ 0.7) bulge-to-total threshold is 
adopted for selecting elliptical galaxies. The distribution of 
stellar masses for elliptical galaxies preserves its bimodal be- 
haviour when considering all galaxies (i.e. is not limited to 
cluster ellipticals). 

• Since ellipticals are the dominant population among 
massive cluster members, one finds that these galaxies have 
been accreted, on average, onto the cluster when residing 
in relatively massive structures (more massive than those of 
the overall cluster galaxy population). This creates a corre- 
lation between the observed fraction of ellipticals and the 
accretion history of the halo that is, however, not strong. 

In the framework of our models, the bimodal distribu- 
tion of elliptical stellar masses is not due to a more promi- 
nent role played by disk instab ility and/or disk regr owth 
for intermediate mass galaxies l|De Lucia et alj l201ll ) . We 
argue that this bimodality results simply from the convolu- 
tion between the strongly decreasing number of galaxies and 
the increasing probability of experiencing a major merger 
event in the past, for increasing ga laxy mass (|De Lucia et al.l 
2006; I Wang fc Kauffmannl 120081 ) . For galaxies with mass 
~ 10 9 M©, the probability of having experienced a major 
merger is not large, but there are many low-mass galaxies. 
On the other hand, almost all galaxies with mass > 10 11 M 
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Figure 3. Mass distribution of haloes accreted on the main branch of each cluster in our sample. Thin dashed lines (blue in the online 
edition of the Journal) are for haloes with an elliptical fraction that is lower than the 10th percentile of the distribution, while thick 
solid lines (red in the online edition of the Journal) correspond to the haloes whose elliptical fraction is larger than the 90th percentile 
of the distribution. The differential distributions in the top panels have been weighted by mass, in order to remove the dominant mass 
dependence and emphasize the differences between the two samples. 



have experienced at least one major merger during their life- 
time so, although the number of massive galaxies is low, this 
mass bin is dominated by elliptical galaxies. 

We argue that clusters that host larger fraction of ellip- 
ticals have a lower contribution from diffuse accretion than 
clusters with lower elliptical fractions (i.e. they accrete more 
haloes, over the entire mass range probed by our simula- 
tions). In addition, we find that the correlation between the 
observed fraction of elliptical galaxies and the accretion his- 
tory of the halo can be weakened in the case of short merger 
times. This would reduce the number of elliptical cluster 
members by having them accreted onto the central clus- 
ter galaxies or merged with other cluster members. In this 
framework, elliptical satellites have been formed before their 
accretion onto the cluster. The measured fraction of ellipti- 
cals is determined by the balance between the disappearance 
of ellipticals due to accretion and mergers, and the recent 
accretion of relatively massive structures (that would likely 
host an elliptical central galaxy). A better 'memory' of the 
accretion history is preserved when merger times are longer. 
In this case, a stronger correlation between the morpholog- 



ical mix of cluster populations and the dynamical status of 
the cluster is expected. 
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